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Abstract 

In  stochastic  acoustic  emission,  both  theory  and  experiments  suggest  that 
the  power  of  the  acoustic  emission  signal  is  proportional  to  the  source  en¬ 
ergy  .  Hence,  inference  about  the  power  is  equivalent  to  inference  about  the 
source  energy  except  for  a  constant  multiple.  In  this  regard,  the  connection 
between  peaks  exceeding  a  fixed  level  and  the  power  in  random  acoustic  emis¬ 
sion  waves  is  explored  when  the  source  energy  is  an  impulse  of  short  duration. 
Under  certain  conditions,  the  peak  distribution  is  sensitive  to  power  changes, 
determines  it  and  is  determined  by  it.  The  maximum  likelihood  estimator 
of  the  power  from  a  random  sample  of  peaks-f/ie  peak  estimatorAs  more  ef¬ 
ficient  than  the  maximum  likelihood  estimator-averaye  sum  of  squares-ixom 
a  random  sample  of  the  same  size  of  signal  values.  When  evaluated  from 
nonrandom  samples,  indications  are  that  the  peak  estimator  may  still  have 
a  relatively  small  mean  square  error.  A  real  data  example  indicates  that  the 
left-truncated  Rayleigh  probability  distribution  may  serve  as  an  adequate 
model  for  high  peaks. 


Abbreviated  Title;  “Acoustic  Emission” 

Key  words  and  phrases;  Energy,  peak  distribution,  Upcrossings,  Maxi¬ 
mum  likelihood,  Gaussian  process,  spectrum. 
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1  Introduction 

In  the  phenomenon  of  acoustic  emission,  a  burst  of  elastic  waves  is  radiated 
into  a  body  from  a  source  undergoing  a  rapid  change  in  stress.  The  acoustic 
emission  problem  consists  of  locating  the  source  and  characterizing  it.  One 
of  the  most  fundamental  parameters  for  source  characterization  is  the  power 
emitted,  a  quantity  proportional  to  the  product  of  the  source  volume  and  ve¬ 
locity.  Experimentally  it  has  been  found  that  the  detectability  of  such  sources 
is  given  by  the  mechanical  power  density  at  the  sensor,  i.e.  the  elastic  power 
of  the  source  divided  by  the  volume  of  the  medium.  Thus  acoustic  emission 
sources  are  undetectable  if  they  are  too  small,  move  too  slowly,  or  are  diluted 
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in  too  large  a  body.  As  in  electromagnetic  theory,  the  dilution  process  oc¬ 
curs  as  the  source  elastic  power  flows  in  the  direction  of  its  acoustic  Poynting 
vector  toward  neighboring  regions  of  lower  power  density.  This  occurs  along 
with  repeated  wave  reflections  from  the  outer  boundaries  in  the  process  of 
reverberation.  During  this  period,  the  waves  tend  to  become  randomized  and 
acoustic  emission  can  be  treated  as  a  stochastic  process.  In  this  paper  we 
examine  source  characterization  for  such  a  process,  beginning  with  a  physical 
model  of  such  a  vibratory  system  adopted  from  seismology.  We  also  examine 
the  accuracy  of  statistical  tools  for  source  energy  characterization. 

Some  important  contributions  have  been  previously  made  to  this  problem. 
Although  Lyon  (1964)  in  his  theory  of  statistical  energy  analysis  adopted  sta¬ 
tistical  methods  to  the  problem  of  energy  transfer  between  components  of  a 
vibrating  system,  some  of  his  techniques  may  be  applicable  to  wave  motion 
within  a  single  body.  For  example  it  has  been  suggested  that  during  surface 
mode  conversion,  energy  is  exchanged  between  the  different  mode  types  ac¬ 
cording  to  statistical  energy  analysis.  Weaver  (1988)  examined  diffuse  wave 
fields.  These  are  fields  of  approximately  constant  energy  which  exist  for 
certain  ranges  of  signal  duration,  specimen  geometry  and  frequency.  Exper¬ 
imentally,  others  (Clough(1987,  1992),  Muravin  et  al.  (1991))  have  found 
that  even  if  frequencies  outside  of  the  diffuse  field  range  are  included  in  the 
measurement  (so  that  the  total  acoustic  energy  is  measured),  the  energy  in 
the  sensor  signal  is  proportional  to  the  source  energy.  This  being  so,  any  es¬ 
timate  of  the  power  from  the  sensor  voltage  is  also  an  estimate  of  the  source 
energy  up  to  a  constant  multiple.  If  the  proportionality  constant  is  known  by 
calibration,  this  linearity  permits  direct  measurements  of  acoustic  emission 
source  energy,  an  example  being  the  energy  released  by  an  advancing  crack. 

One  of  the  principal  features  examined  in  this  study  is  the  relationship 
for  such  systems  between  power  measured  directly  from  the  signal  versus 
that  obtained  from  the  distribution  of  its  larger  peaks  (truncated  signal). 
The  relationship  between  the  power  and  peak  magnitude  and  its  probability 
distribution  truncated  to  the  /e/t-conditional  on  peak  magnitude  greater  than 
some  threshold-is  the  main  topic  of  the  present  paper.  To  gain  insight  as 
well  as  some  generality,  the  underlying  process  is  assumed  to  be  a  monotone 
transformation  of  a  Gaussian  stationary  process.  The  Gaussian  case  is  then 
a  special  case  only. 

Specifically,  in  this  paper  it  is  shown  that: 
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1)  The  peak  distribution,  conditional  on  peak  size  greater  than  a  fixed 
positive  threshold  L,  depends  on  the  power.  Changes  in  the  conditional 
peak  distribution  correspond  to  changes  in  power. 


2)  Under  certain  conditions,  the  maximum  likelihood  estimator  of  the  power 
from  a  random  sample  Xi,...,Xn  of  peak  heights  greater  than  a  rela¬ 
tively  large  L  >  0, 


£"=i  Xf  -  nL-‘ 

2n 

is  more  efficient  than  the  more  direct  estimator  obtained  from  the  sum 
of  squares  of  a  random  sample  Zi, ...,  of  signal  values. 


(T 
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3)  When  and  are  computed  from  a  time  series  (i.e.  not  a  random 

sample),  computer  simulations  indicate  that  there  are  levels  L  for  which 
tends  to  have  a  smaller  mean  square  error  (MSE)  than 

4)  The  conditional  peak  distribution  of  any  process  obtained  via  a  mono¬ 

tone  transformation  of  a  stationary  Gaussian  process,  constitutes  a  one 
parameter  exponential  family  which  generalizes  the  Rayleigh  distribu¬ 
tion. 

5)  Power  is  a  feature  useful  in  detecting  changes  in  the  source  or  in  the 

medium. 


2  Preliminaries 

The  existence  of  a  basic  linear  relationship  between  source  energy  and  average 
power  in  the  sensor  signal  was  discussed  above.  This  is  what  triggered  the 
present  investigation. 

2.1  An  Energy  Model 

Experimental  results  show  that  when  an  energy  impulse  of  short  duration  r, 
say,  operates  on  a  metal  plate,  the  expected  energy  density  S{t)  in  the  plate 


3 


rises  sharply  throughout  the  period  (0,  r],  and  then  dissipates  gradually  for 
t  >  T.  The  latter  period  is  referred  to  as  ringdown  or  reverberation.  A  basic 
physical  model  which  describes  this  behavior  runs  as  follows  (e.g.  see  Wood 
(1966),  pp.  536-8).  For  simplicity  we  assume  that  the  energy  impulse  is  of 
constant  intensity  Vs  in  (0,  r].  Let  m  >  0,  be  an  absorption  constant  which 
depends  on  the  geometry  and  material  of  the  specimen.  Then,  as  in  Clough 
(1987),  for  0  <  <  <  T  and  with  £(0)  =  0,  the  time  derivative  is 

£'{t)  —  Vam£(t) 

so  that  the  expected  energy  density  increases  as 

(1)  £{t)  =  —(1  —  exp(— mt)),  0  <  <  <  r 

TtX 

Since  the  source  ceases  operating  for  t  >  t  we  have 

£'{t)  =  —m£{t) 

and  because  now  £(t)  =  Vs{l  —  exp(— mr))/m,  we  obtain  for  the  ringdown 
period  the  solution 

V 

(2)  £(t)  =  —(1  —  exp(— mT))exp(— m(t  —  t)),  t  >  t 

m 

Now  fixing  arbitrary  time  points  ti,t2  >  r,  we  observe  that  the  power  over 
[^1,^2]  is  given  by 

f  *  £{t)dt  =  I  ?"^’^\(exp(-mti)  -  exp(-mf2))|  Vs 
t2  — 1\  Jti  ( ni(i2  —  *1)  J 

where  C(m,T)  is  a  constant  depending  only  on  m  and  r.  It  follows  that 
the  source  energy  is  proportional  to  the  power  over  [^1,^2]  for  any  ti  <  ^2- 
Moreover,  we  can  choose  a  sufficiently  small  interval  [^1,^2]  for  which  the 
mechanical  vibration  is  approximately  stationary. 

Note  that  we  have  not  yet  touched  on  the  observed  electrical  signal.  This 
is  done  next. 
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2.2  A  Statistical  Vibration  Model 


The  mechanical  vibration  at  a  location  on  the  plate  is  closely  related  to  the 
vibration  at  other  locations.  In  other  words,  a  time  series  representing  the 
vibration  at  a  certain  location  as  a  function  of  time  depends  on  similar  series 
observed  at  other  locations.  If  we  assume  that  the  different  time  series-their 
number  can  be  very  large-are  related  to  each  other  linearly,  their  mutual 
interdependence  may  be  approximated  by  a  first  order  linear  system  of  dif¬ 
ferential  equations.  Thus,  from  a  purely  mathematical  point  of  view,  an 
idealized  model  which  suits  our  purpose  has  the  form 

(3)  =  Ay(t)  +  h6{t) 

with  the  initial  condition  y(0— )  =  0.  Here  the  displacement  y(t)  is  an  A;  x  1 
column  vector  and  A  is  a  real  square  matrix  of  degree  k  x  k.  Bolt  and 
Brillinger  (1979)  and  Brillinger  (1987)  note  that  (  3)  serves  as  an  adequate 
model  for  a  great  variety  of  vibratory  motions. 

The  general  solution  of  (  3)  is  (Boyce  and  DiPrima  (1986),  Ch.  7,  Cod- 
dington  and  Levinson  (1955),  p.  78), 

(4)  y{t)  =  +  h  f  e^^^~‘H{s)ds  =  t  >  0 

no¬ 

where 

oo  .k 

e*‘  =  H-EA‘|; 

Jfc=l 

Assume  the  matrix  A  has  n  distinct  complex  eigenvalues  \j.  Then  the 
corresponding  eigenvectors  (column)  Xj  are  complex  and  are  linearly  inde¬ 
pendent.  If  we  define  the  matrix 


T  =  (Xi,X2,...,X„) 

then  it  is  nonsingular.  Now,' 

AT  A(xi ,  X2, . . . ,  Xn)  ( Axj,  Ax2,  ...,  Axn)  —  (A^xi,  A2X2, ...,  AnX„) 

=  (xi,X2,...,x„)diag(Ai,A2,...,An)  =  TD 

so  that  A  is  diagonalizable, 

T-*AT  =  D 
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or 

A  =  TDT-^ 

Substituting  this  in  we  obtain 

oo  j^k 

gAi  =  T{I  +  X]  =  Te^'T"^ 

A:=l 

But 

e°‘  =  diag(e'^^*,e^^*,...,e^"*) 
and  our  solution  (  4)  becomes 

(5)  y{t)  =  =  Tdiag(e^'',  ^ 

j=i 

where  aj  is  the  jth  component  of  the  vector  T~^b. 

Now  write 

Aj  =  <Tj  -j-  iu>j,  Xj  =  a.j  -b  ibj 
where  the  ringdown  effect  suggests  that  cTj  >  0.  Then 

e’^-’*Xj  =  cosojjt  —  bj  sinwji}  +  ie~'^^*{hj  coswjt  +  aj  sinwjt} 

Multiplying  this  by  (complex)  aj,  it  follows  that  each  component  of  y{t),  say 
y{t),  has  a  real  solution  of  the  form 

(6)  y(t)  =  ^  cos(a;j^  +  <fj) 

j 

where  — <Tj  =  3JAj,  Wj  = 

The  solution  (  6)  indicates  that  the  time  series  of  mechanical  vibration 
at  a  point  resulting  from  acoustic  emission  may  be  modeled  for  a  finite  K  as 

K 

(7)  y(t)  =  cos{ujt  +  (fj)  +  e{t) 

j=i 

where  {e(f)}  is  noise  (see  Brillinger  (1987)).  Evidently  (  7)  is  conducive  to 
explaining  the  ringdown  effect.  When  the  amplitudes  and  phases  are  uncor¬ 
related,  the  signal  to  noise  ratio  is  relatively  high,  and  a  dominant  frequency 
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exists  (  see  Figures  3  and  4),  the  expected  energy  in  (  7)  follows  approxi¬ 
mately  the  pattern  of  the  energy  model  (  2).  Estimation  of  the  parameters 
in  (  7)  is  discussed  in  Bolt  and  Brillinger  (1979). 

A  transducer  which  transforms  the  mechanical  vibration  into  the  observed 
electrical  signal,  acts  as  a  linear  filter  on  {2/(<)}.  Carrying  our  ideas  a  bit 
further,  if  the  aj  in  (  7)  are  sufficiently  small,  then  the  observed  electrical 
signal  over  a  short  interval,  denoted  by  {Z{t)},  may  be  taken  as  a  mixed 
spectrum  stationary  process  with  mean  0  and  variance  of  the  form  (see 
also  Rice  (1944,1945)) 

K' 

(8)  Z{t)  =  Y,  ^3  C.Os{u3jt  -I-  (f>j)  +  C(t) 

j=i 

with  the  damping  factors  removed.  Thus,  by  the  mean  value  theorem  and 
in  light  of  the  preceding  subsection,  the  variance  cr^-that  is,  the  power  of 
{Z{t)}-is  proportional  to  Va-  The  experimental  results  in  Clough  (1987,1992) 
support  this  conclusion. 

2.2.1  Experimental  Verification 

To  see  whether  (  8)  is  a  reasonable  model  for  acoustic  emission,  we  have 
generated  acoustic  emission  signals  in  the  laboratory  and  computed  their 
FFT’s  and  autocorrelations.  The  data  and  the  method  of  generation  are 
described  next. 

Acoustic  emission  signals  were  produced  by  releasing  a  3.175  mm  diameter 
steel  ball  bearing  from  a  fixed  height  of  40.5  cm  above  a  plate  of  2519-T87 
hardened  aluminum  alloy.  The  plate  had  dimensions  of  2.22  x  30.48  x  91.44 
cm,  supported  at  each  corner  by  spheres.  For  some  of  the  experiments  a  block 
of  2024-T351  hardened  aluminum  alloy  with  dimensions  11.43  x  8.25  x  8.25 
cm  was  placed  on  the  plate  center  ,  with  its  large  face  down.  It  was  coupled 
to  the  plate  with  ultrasonic  gel. 

The  locations  of  the  impact  point  of  the  ball  and  the  transducer  were 
along  the  centerline  of  the  plate,  15  cm  from  each  end.  The  acoustic  emission 
transducer  was  a  piezoelectric  type,  coupled  by  an  ultrasonic  gel  to  the  plate 
and  loaded  with  a  0.6  N  steel  weight  to  maintain  reproducibility.  The  cali¬ 
bration  supplied  by  the  manufacturer  showed  a  dominant  resonant  response 
for  the  transducer  at  140  kHz.  The  transducer  was  coupled  by  a  coaxial  ca¬ 
ble  directly  into  the  high-impedance  input  of  a  12  bit  digital  recorder  with  a 
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preset  threshold  voltage  for  triggering  the  recording.  Each  recording  of  4096 
points  was  made  with  a  200  ns  sampling  time  interval,  and  the  recorder  per¬ 
mitted  recording  points  before  triggering  occurred.  The  voltages  recorded 
typically  had  a  maximum  range  of  ±  2  volts,  corresponding  to  an  integer 
range  of  ±2048. 

Figures  1  and  2  show  segments  of  typical  electrical  signals  obtained  by 
dropping  the  small  ball  bearing  on  the  metal  plate  and  also  on  the  same 
plate  plus  a  metal  block  shortly  before  a  pronounced  ringdown  takes  place. 
The  qualitative  appearance  of  the  records  make  the  stationary  assumption 
over  a  short  interval  more  credible.  The  corresponding  averaged  (over  several 
realizations)  FFT’s  in  Figures  3  and  4  show  that  the  data  contain  strong 
sinusoidal  components  in  close  agreement  with  the  model  (  8).  Figures  5 
and  6  show  the  corresponding  estimated  autocorrelation  functions  which  do 
not  seem  to  decay-consistent  with  the  presence  of  line  spectra.  This  means 
that  the  statistical  dependence  between  well  separated  observations  does  not 
vanish.  The  estimation  of  power  under  this  condition  can  be  quite  tricky 
(e.g.  see  Yaglom  (1987),  pp.  226-7,  Koopmans  (1974)  p.  60). 

The  relationship  between  the  source  and  measured  acoustic  emission  en¬ 
ergies  produced  by  the  steel  ball  bearing  dropped  from  various  increasing 
heights  onto  the  aluminum  alloy  plate  is  shown  in  Figure  7.  With  a  sample 
correlation  coefficient  of  .974,  the  observed  relationship  is  linear  to  a  large 
degree,  in  line  with  our  energy  model  which  predicts  a  highly  linear  rela¬ 
tionship.  A  least  squares  fit  of  the  model  Y  =  aX^  to  the  data  gives  the 
estimates  a  =  0.8979  and  b  =  1.27165.  The  slight  curvature  is  associated 
with  an  energy  dissipation  mechanism  in  the  plate  at  the  impact  point,  since 
the  measured  rebound  height  is  much  less  than  that  predicted  by  the  Hertzian 
elastic  theory  of  impact.  This  theory  has  been  verified  experimentally  by  Ra¬ 
man  (1920)  as  well  as  by  others.  Two  such  mechanisms,  which  would  also 
produce  the  observed  nonlinear  response,  are  plastic  deformation  in  the  plate 
or  fracture  of  the  aluminum  oxide  at  the  plate  surface. 


3  Power  Estimation  Prom  Peak  Distribution 

The  relevance  of  the  probability  distribution  of  peak  magnitude  (measured 
from  level  0)  in  reliability  studies  of  mechanical  systems  has  been  known  for  a 
long  time  (e.g.  see  Lin  (1967)  Ch.  9,  Yang  (1986)  Ch.  8).  Regarding  acoustic 
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emission,  Ono  (1976)  studied  the  peak  distribution  of  continuous  and  burst- 
type  acoustic  emission  signals.  An  argument  is  made  there  that  the  Rayleigh 
distribution  serves  as  an  adequate  model  for  the  peak  magnitude  observed 
in  continuous-type  acoustic  emission  regardless  of  bandwidth.  The  work  in 
these  references  is  based  to  a  large  degree  on  the  work  of  Rice  (1944,  1945) 
and  the  Gaussian  assumption.  Builo  and  Tripalin  (1982)  have  demonstrated 
experimentally  that  the  exponential  distribution  is  an  adequate  probability 
distribution  model  for  the  intervals  between  high  peaks  of  acoustic  emission 
observed  in  low-carbon  steel. 

What  is  of  interest  to  us  is  to  further  explore  the  connection  between 
power  and  peak  distribution  conditional  on  peaks  higher  than  a  fixed  pos¬ 
itive  threshold.  From  the  preceding  discussion,  this  enables  statistical  in¬ 
ference  about  the  relative  (i.e.  up  to  a  constant)  source  energy  impulse.  In 
what  follows,  the  observed  electrical  signal  is  assumed  to  be  a  coordinatewise 
monotone  transformation  of  a  Gaussian  process. 

3.1  A  Generalization  of  the  Rayleigh  Distribution 

We  shall  derive  the  conditional  peak  distribution  in  monotone  transforma¬ 
tions  of  Gaussian  processes  to  find  the  sensitivity  of  the  distribution  to 
changes  in  power.  Observe  that  a  monotone  transformation  of  a  Gaussian 
process  is  not  necessarily  Gaussian,  a  fact  which  makes  our  results  somewhat 
more  general.  In  what  follows  we  always  assume  that  the  relevant  spectral 
moments  exist. 

Let  {Z{t)}j  — oo  <  t  <  oo,  E[Z(<)]  =  0,  Var[Z(f)]  =  cr^,  be  a  stationary 
Gaussian  process  with  autocorrelation  Pz{t),  assumed  twice  differentiable  at 
f  =  0.  The  process  in  (  8)  fulfills  this  assumption  when  the  amplitudes  are 
independent  (positive)  Rayleigh  random  variables,  the  phases  are  indepen¬ 
dent  random  variables  uniformly  distributed  in  [0, 27r]  independently  of  the 
amplitudes,  and  the  noise  is  a  zero  mean  continuous  spectrum  stationary 
Gaussian  process  independent  of  the  amplitudes  and  the  phases.  Consider 
the  memoryless  transformation  of  Z(t), 

Y{t)  =  cp{Z(t)) 

where  (p(x)  is  strictly  monotone  increasing  and  continuously  differentiable. 
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Observe  that  <7^  is  the  power  of  {Z{t)},  while  that  of  can  be 

approximated  roughly  from  a  Taylor’s  series  expansion  to  one  term, 

Var[y(<)]  (^'(0))^cr^ 

A  much  more  refined  approximation  obtained  at  the  price  of  higher  deriva¬ 
tives  is  obtained  from  Houdre  and  Kagan  (1992).  Let  X  ~  A/’(0,  <t^),  and 
(p{x)  be  a  complex-valued  2n  times  differentiable  function  such  that 

<  oo,  A:  =  0, 1,  •  •  • ,  2n,  n>l 

Then  the  Houdre-Kagan  inequality  states  that 

2n  /  lU+l 

E  <  VarMX)]  <  E 

k=l  k=l 

By  some  integration,  the  joint  probability  density  of  {Y (t),  g{u,  v), 

is  given  by  (Orsingher  (1979)), 

la\  «l.,  (->;(°))M>-»(«))p]} 

u  >  </?(— oo),  — oo  <  V  <  oo.  Following  Rice  (1944,  1945),  the  expected  rate- 
per  unit  time-of  upcrossings  of  level  u  for  {y(f)},  denoted  by  Ny{u),  is  then 
the  integral  (Orsingher  (1979)) 

(10)  Ny{u)  =  J^  vg{u,v)dv  =  ^^^-^"(0)  exp  |^[(/p-^(u)]2| 

where  u  >  (p(—oo).  Related  formulas  both  in  discrete  and  continuous  time 
are  discussed  in  Barnett  and  Kedem  (1991),  and  for  a  general  discussion  see 
Kedem  (1994),  Ch.  4. 

By  a  peak  we  shall  mean  a  local  positive  maximum  observed  in  {y(i)}. 
Note  that,  in  a  given  time  interval,  the  number  of  peaks  located  above  a 
sufficiently  high  threshold  level  L{>  0)  is  equal  to  the  number  of  upcrossings 
of  L  under  the  condition  that  there  are  no  troughs  above  L.  This  simple  fact 
enables  us  to  derive  the  probability  distribution  of  large  peaks  from  (10). 

Denote  by  the  random  variable  X  the  magnitude  of  a  peak  observed 
in  {y(^)}.  From  now  on  we  assume  that  L  is  a  sufficiently  large  positive 
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constant  so  that  the  probability  of  a  trough  above  L  is  negligible.  Then  from 

(10) 


P{x  <  X  <  X  +  dx\X  >  L)  = 


Ny(x)  —  Ny{x  +  dx) 


—Ny(x)dx 

N.(l) 


from  which  we  obtain  the  conditional  density  of  X  given  that  X  >  L  as 

(11)  p(s:.  c)  =  exp  ^^l(^-'(x)f  -  (.^-‘(i))^]}  ,  x  >  £ 

We  can  see  that  (  11)  defines  a  one  parameter  exponential  family  in  a  for 
any  given  monotone  <f{x)  and  fixed  L  >  0,  the  sufficient  statistic  being 
{(p~^(X)y.  For  (p(x)  =  X  and  T  =  0,  p{x,cr)  reduces  to  the  well  known 
Rayleigh  probability  density  function  (pdf), 

(12)  p((c,<7)  =  ^expj^j,  a;  >  0 


Evidently,  (11)  provides  a  generalization  of  the  Rayleigh  distribution. 

To  illustrate  how  changes  in  the  parameter  cr  change  the  general  pdf 
p{x,(j)  in  (  11),  we  have  plotted  p{x,(t)  in  Figures  8  to  13  using  various 
monotone  transformations  and  different  levels  L.  The  transformations  are 
'p{x)  —  a:,  exp(a;),  1/(1 +  exp(— x)),arctan(x).  We  can  see  from  the  figures 
that  the  conditional  peak  distribution  is  fairly  sensitive  to  changes  in  cr-a 
fact  which  might  suggest  a  or  as  a  classification  feature  in  characteriz¬ 
ing  acoustic  emission  types.  Some  evidence  supporting  this  is  derived  from 
the  acoustic  emission  data  (ball  dropping)  described  earlier.  Thus  from  9 
independent  realizations  of  “plate”  only  and  then  from  10  independent  real¬ 
izations  of  “plate-fblock”  we  obtain  the  sample  variances  (normalized)  as  de¬ 
picted  in  Table  1.  The  acoustic  emission  variance  (power)  for  “plate+block” 
tends  to  be  appreciably  lower.  Crude  95%  confidence  intervals  for  “plate” 
and  “plate-fblock”  power  are  (0.608,0.913)  and  (0.317,0.523),  respectively. 
The  intervals  do  not  intersect. 

It  is  useful  to  note  that  the  mean  and  variance  of  the  sufficient  statistic 
are  readily  available  once  p{x,a)  is  reparametrized  using  p  = 
Then  we  obtain  the  one  parameter  exponential  family  in  natural 
form  (Bickel  and  Doksum  (1977),  pp.  70-71), 

p{x,cr)  =  exp|77((,c>-^(x))2  -|-log(-277)  -7/(v?"^(T))2  +  log[^;^^^^| 
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With 

do{ri)  =  log(-277)  -  ^{L))^ 

we  have 

(13) 

EJ(^-^(X))^]  =  -dM  =  2a^  +  {^-\L)) 

and 

(14) 

Since  inference  in  exponential  families  is  based  on  the  sufficient  statistics, 
the  moment  expressions  (  13)  and  (  14)  are  very  helpful. 
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Table  1.  Estimated  power/lO^  for  acoustic  emission  records  obtained  from 
“plate”  and  “plate+block” . 


Plate 

Plate-fB 

0.71528 

0.45827 

0.73640 

0.66554 

0.40252 

0.75503 

1.10723 

0.28555 

0.94713 

0.37602 

0.68071 

0.39467 

0.76387 

0.29357 

1.01470 

0.24337 

0.47773 

0.37462 

0.35065 

Mean  0.76062  0.41973 

SD  0.23366  0.16623 


3.2  Estimation  of  Power  Prom  Peaks 

By  our  construction,  the  estimation  of  the  power  in  {T(i)}  is  equivalent 
to  the  estimation  of  when  'p{x)  is  known-an  assumption  we  have  made 
implicitly  all  along.  Thus,  the  estimation  can  be  carried  out  in  several  ways 
of  which  we  shall  single  out  estimation  based  on  actual  values  and  estimation 
based  on  peaks  larger  than  L{>  0).  It  is  interesting  to  compare  the  estimates 
constructed  in  these  two  ways.  As  we  show  below,  there  can  be  an  advantage 
in  precision,  and  possibly  in  speed,  in  estimating  power  from  peak  heights 
rather  than  from  sequential  voltages. 

As  noted  above,  acoustic  emission  data  may  contain  discrete  spectral  com¬ 
ponents  so  that  the  statistical  dependence  in  the  data  never  dies  out  when 
observations  are  separated  in  time.  In  the  Gaussian  case  we  have  the  further 
complication  that  the  sample  variance  is  altogether  inconsistent  in  the  pres¬ 
ence  of  line  spectra.  Thus  the  power  estimation  problem,  which  ostensibly 
appears  to  be  rather  simple,  requires  some  caution.  For  this  reason,  we  shall 
first  assume  that  one  has  at  hand  random  samples  of  independent  and  iden¬ 
tically  distributed  random  observations  from  the  process  values  and  from  its 
peaks.  In  source  characterization  and  reliability  studies,  the  scientist  may 
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stimulate  a  structure  to  generate  many  independent  acoustic  emission  records 
by  impact  or  by  using  a  laser  source  at  a  very  low  cost.  Hence,  the  assumption 
that  random  samples-extracted  from  acoustic  emission  records-are  available 
has  a  practical  base.  The  random  sample  approach  enables  maximum  likeli¬ 
hood  considerations  which  in  many  cases  result  in  efficient  estimators.  This 
also  may  suggest  a  general  form-at  times  quite  unexpected-for  the  estimator 
regardless  of  statistical  independence. 

So,  for  a  relatively  large  L{>  0),  let  Xi,X2,...,X„  be  a  random  sample 
of  peaks  greater  than  L  from  {y(f)}.  Then  from  (  11)  the  sufficient  statistic 
for  <T  is 

T  = 

i=l 

and  from  (  13), (  14), 

E^[T]  =  2na^  +  n{(p~^{Lyf 

and,  regardless  of  ^{x), 

V&Xa  [T]  =  4ncr^ 

The  maximum  likelihood  estimator  (MLE)  of  is  obtained  from  the  equa¬ 
tion  (Bickel  and  Doksum  (1977),  p.  102) 

2n(T^  +  n{ip-^{L)f  = 

t=i 

which  gives 

2n 

It  is  easily  seen  that  cr^  is  unbiased  and 

(16)  Var,[<;2]  ^  ^ 

n 

Moreover, 

—  cr^)  A/’(0,  (7^) 

Suppose  on  the  other  hand  we  have  a  random  sample  of  values  Yi^Yz,  ...,Yn 
from  {y(t)}.  Then  the  maximum  likelihood  estimator  of  cr^  becomes 


(17) 


n 
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which  is  again  unbiased,  but 


(18)  Var^[<T2]  = 

n 

which  is  twice  as  large  as  the  MLE  from  (16).  We  therefore  have 

Theorem  3.1  The  MLE  <7^  from  a  random  sample  of  size  n  of  peaks  exceed¬ 
ing  a  positive  threshold  L  is  twice  as  efficient  as  the  MLE  cr'^  from  a  random 
sample  of  size  n  of  actual  process  values. 

It  follows  that  we  must  have  twice  as  many  actual  independent  process  values 
than  independent  peaks  to  achieve  the  same  precision. 

We  need  to  clarify  the  last  result.  Recall  that  we  are  dealing  with  an  ideal 
situation  where  potential  difficulties  stemming  from  discrete  spectral  compo¬ 
nents  are  bypassed  by  sampling  from  independent  realizations.  Then  under 
some  conditions  the  power  can  be  estimated  efficiently  from  large  peaks. 
However,  the  presence  of  line  spectra  is  not  necessarily  a  hindrance  and  we 
can  still  use  the  common  sense  estimator  which  is  the  average  of  the  sum  of 
squares  from  a  single  realization-a;SSuming  zero  mean.  For  example,  when 
in  (  8)  the  amplitudes  are  fixed  (not  random),  the  phases  are  uniformly  and 
independently  distributed  over  the  spectral  support  (0, 27r]  (assuming  the 
unit  of  time  to  be  the  sampling  interval),  and  the  noise  component  is  ergodic 
and  independent  of  the  random  phases,  then  the  common  sense  estimator 
is  strongly  consistent.  For  discussions  of  this  point  see  Houdre  and  Kedem 
(1994),  Li  and  Kedem  (1993),  Kedem  and  Slud  (1994),  and  Kedem  (1994), 
Ch.  6.  A  clarification  of  this  technical  point  is  given  in  the  enclosed  ap¬ 
pendix.  Whether  in  this  case-the  fixed  amplitudes  case-it  is  more  efficient 
to  estimate  the  power  from  peaks  is  not  at  all  clear  (although  our  simulation 
study  below  clarifies  this  point  to  some  degree),  but  one  thing  is  certain, 
there  are  many  more  sample  values  than  qualified  large  peaks  exceeding  L 
in  a  single  realization,  a  fact  which  may  make  the  use  of  peaks  somewhat 
dubious.  Still,  the  import  of  our  finding  is  that  it  points  to  the  connection 
between  power  and  peak  values,  and  to  the  advantage  of  peaks  under  some 
conditions.  Furthermore,  as  we  shall  very  soon  see,  the  estimator  (  15)  can 
be  used  advantageously  in  general  regardless  of  statistical  independence. 
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4  The  Case  of  Truncated  Rayleigh 


We  now  specialize  the  preceding  general  discussion  to  the  case  of  most  prac¬ 
tical  interest  when  no  transformation  takes  place  at  all-that  is  <p{x)  =  a;-and 
the  power  cr^  is  estimated  from  large  peaks  in  the  Gaussian  process  {Z{t)}. 
That  is,  we  assume  that  the  observed  electrical  signal  is  stationary  and  Gaus¬ 
sian  with  mean  zero.  Then  (10)  becomes 

(19)  Ny(^)  = 

When  the  measurement  level  u  is  sufficiently  high  (  19)  is  also  the  expected 
rate  of  peaks  above  level  u  (see  Rice  (1945),  p.  75)  and  the  conditional 
density  (  11)  of  peak  magnitude  X  given  that  X  >  L  fox  L  >  0  becomes 

(20)  p{x,(t)  =  ,x  >  L 


which  is  a  (left-)  truncated  Rayleigh  pdf.  For  T  =  0  no  truncation  occurs  and 
(  20)  reduces  to  a  bona  fide  Rayleigh  probability  density.  The  MLE  from  a 
random  sample  Xi X^  of  peaks  greater  than  L  now  has  the  form 

T.U 

2n 

This  is  the  peak  estimator  ixoxa  which  an  approximate  95%  confidence  interval 
for  the  power  is  given  by 


(22) 


<T^  =  <T^  ±  1.96 


^/n 


When  instead  of  peaks  we  have  a  random  sample  of  process  values  Zi, ..., 
the  MLE  is 

y2 

(23)  0^2  3=  ^*=1  ‘• 


n 


In  this  case  the  approximate  95%  confidence  interval  is  longer. 


(24) 


CT 


=  ^2  ±  1.90 


Since  cr^  is  proportional  to  Vs,  from  a  single  calibration  which  determines 
the  proportionality  constant  we  can  get-using  either  (  22)  or  (  24)-a  confi¬ 
dence  interval  for  Va- 
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4.1  Plate  Data  Example:  Large  Peaks  Are  Rayleigh 

To  see  whether  (  20)  is  a  reasonable  model  for  acoustic  emission  signals,  the 
truncated  density  was  fitted  to  the  data  from  our  plate  experiment  (plate 
only)  described  in  Section  2.2.  One  may  argue  that  high  peaks  are  roughly 
nearly  independent  and  we  have  used  all  the  peaks  greater  than  L  =  400 
from  9  independent  electrical  signals.  The  results  of  a  chi-square  goodness  of 
fit  test  are  given  in  Table  2.  The  results  are  insignificant  at  level  .05,  that  is, 
the  hypothesis  that  the  peak  data  follow  a  truncated  Rayleigh  with  L  =  400 
is  accepted  at  the  .05  level. 

On  the  other  hand,  lowering  L  may  give  significance.  Thus,  with  L  =  300 
and  an  additional  category,  300  <  X  <  400,  we  have  a  large  observed 
value  as  indicated  in  Table  3. 

Table  2.  Observed  versus  expected  frequencies  under  the  truncated  Rayleigh 
hypothesis  with  L  =  400.  The  hypothesis  is  accepted  at  level  .05. 


Category 

Observed 

Expected 

400  <  X  <  500 

40 

28.93 

500  <  A  <  600 

23 

24.87 

600  <  X  <  700 

12 

19.26 

700  <  X  <  800 

13 

13.58 

X  >  800 

18 

19.37 

=  1.41  X  10^ 
=  7.2 
d.f=3 


=  7.8 
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Table  3.  Observed  versus  expected  frequencies  under  the  truncated  Rayleigh 
hypothesis  with  L  =  300.  The  hypothesis  is  rejected  at  level  .05. 


Category 

Observed 

Expected 

300  <  X  <  400 

83 

53.44 

400  <  a:  <  500 

40 

47.14 

500  <  a:  <  600 

23 

35.97 

600  <  a:  <  700 

12 

24.16 

700  <  X  <  800 

13 

14.41 

X  >  800 

18 

13.88 

(P'  =  1.05  X  10® 
=  29.6 
d.f=4 


xL  =  9-48 

We  conclude  that  the  truncated  Rayleigh  model  (  20)  with  a  “large”  L 
is  a  reasonable  model  for  our  acoustic  emission  signals.  Note  that  assuming 
a  truncated  Rayleigh  model  for  the  peaks  is  weaker  than  the  full-fledged 
Gaussian  assumption. 

As  for  the  peak  estimator  (  21),  a  sample  of  30  peak  values  above  L  =  750 
from  the  plate  data  gives  cf^  —  1.01790  x  10®,  which  is  in  agreement  with 
some  of  the  higher  values  in  the  “Plate”  column  in  Table  1.  The  peaks  larger 
than  L  =  650  from  a  single  ball  dropping  are  992, 656, 768, 704, 672, 672, 800 
and  =  0.77536  X  10®,  a  figure  close  to  the  average  of  the  sample  variances 
in  the  “Plate”  column  in  Table  1.  However,  the  samples  here  are  not  exactly 
random,  and  their  size  is  too  small  for  a  meaningful  comparison  with  the 
estimates  in  Table  1  obtained  from  thousands  of  data  points.  More  definite 
conclusions  may  be  reached  from  simulations  under  controlled  conditions  as 
in  the  next  subsection. 


4.2  Comparison  Between  and  cr^ 

4.2.1  Random  Samples 

Consider  the  stationary  Gaussian  autoregressive  process  of  order  two,  Ai7(2), 


Zt  —  (j)xZt-\  -|-  <l>2Zt-2  +  e* 
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where  t  =  ±1,±2,...,  and  the  Ct  are  independent  A/’(0, 1).  We  have  simu¬ 
lated  time  series  of  length  10,000  from  the  process  for  various  choices  of  </>2. 
From  each  stretch  of  100  points-there  are  100  stretches-we  selected  one  pro¬ 
cess  value  and  one  peak  greater  than  a  fixed  L.  The  samples  obtained  in 
this  way  are  close  to  being  random  samples  since  the  autocorrelation  decays 
exponentially  fast  and  the  process  is  Gaussian. 

For  (</>!,  </>2)  =  (0.7,  —0.5),  the  true  variance  of  the  process  is  cr^  =  1.70455. 
Table  4  shows  the  results  of  15  independent  time  series  of  length  10,000  each 
from  the  process,  each  giving  rise  to  values  of  with  L  =  2,  and  cr^.  The 
average,  sample  variance,  and  mean  square  error  (MSE)  from  15  values  show 
that  the  peak  estimator  outperforms  the  average  sum  of  squares  cr^.  It  is 
interesting  to  note  that  the  MSE  for  is  about  one  half  of  that  of  cr^,  in 
agreement  with  the  theory. 

Table  4.  Values  of(cr‘^,a‘^)  corresponding  to  (^1,^2)  =  (0.7, —0.5),  L  =  2, 
n  =  100.  The  true  power  is  =  1.7045. 

^2  c?2 

1.60930  1.74627 
1.32313  1.69699 
1.48443  1.57510 
2.00502  1.84070 

1.50617  1.61886 
1.74592  1.47999 

1.34961  1.56067 
1.56412  1.73637 
1.53982  1.60481 
1.56688  1.47263 
1.47726  1.77821 
1.43776  1.62170 
1.73009  1.72299 

1.78605  1.98918 
1.57093  1.61198 
Average  1.57977  1.67043 
Variance  0.03147  0.01893 
MSE  0.04704  0.02010 

The  same  experiment  was  repeated  with  15  independent  time  series  of 
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length  20,000  each  from  the  process  with  (^1,(^2)  =  (0.3, —0.3).  Here  each 
stretch  of  200  points  from  a  given  series  gave  us  one  process  value  and  one 
peak  above  L  =  2,  so  that  n  =  100.  Again  we  can  see  from  Table  5  a  marked 
decrease  in  MSE  in  favor  of  the  peak  estimator. 

Table  5.  Values  of  corresponding  to  ((^1,(^2)  =  (0.3,— 0.3),  L  =  2, 

n  =  100.  The  true  power  is  =  1.16071. 


<7^ 

<j2 

0.94680 

1.20669 

0.82511 

1.09422 

1.49400 

1.15819 

1.04945 

1.11955 

1.10327 

1.12631 

1.08123 

1.14332 

1.21339 

0.91157 

0.82539 

1.15223 

1.01798 

1.03584 

1.32814 

1.07454 

0.91231 

1.16006 

1.19762 

1.05037 

1.04856 

1.18575 

1.01346 

1.01346 

1.40548 

1.26507 

Average  1.09748  1.11314 
Variance  0.03945  0.00759 
MSE  0.04345  0.00986 


The  preceding  example  was  repeated  for  (^1,^2)  =  (0.5, —0.5)  but  with 
a  relatively  large  L  =  3.  Here  not  every  stretch  of  200  produced  a  peak 
larger  than  T  =  3,  so  that  we  ended  up  evaluating  cr^  from  less  than  n  =  100 
values.  However,  the  average  sum  of  squares  was  still  evaluated  from 
n  =  100  values.  The  results  of  20  runs  are  given  in  Table  6.  The  table  also 
gives  the  number  of  qualified  peaks  used  in  the  evaluation  of  Once  again, 
the  MSE  of  cr^  is  smaller  than  the  MSE  of  cr^,  despite  of  the  fact  that  the 
peak  estimator  was  obtained  from  smaller  samples. 
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Table  6.  Values  of  {cr^^a'^)  corresponding  to  =  (0.5, —0.5)^  L  —  3. 

For  n  =  100.  The  true  power  is  =  1.5. 


0-2 

<7^  No.  Peaks 

1.77575 

1.55380 

81 

1.64088 

1.57567 

77 

1.81782 

1.46814 

81 

1.81136 

1.56932 

70 

1.71470 

1.41060 

72 

1.26101 

1.43511 

71 

1.35314 

1.31085 

65 

1.57529 

1.47282 

80 

1.75162 

1.35096 

74 

1.86053 

1.40739 

77 

1.33553 

1.56781 

66 

1.48608 

1.32215 

78 

1.19041 

1.25009 

70 

1.41774 

1.43028 

77 

1.31050 

1.55496 

74 

1.77570 

1.50617 

68 

1.20113 

1.57117 

72 

1.68260 

1.56711 

74 

1.46441 

1.44517 

76 

1.82231 

1.61418 

79 

Average 

1.56243 

1.46919 

Variance 

0.05312 

0.01090 

MSE 

0.05702 

0.01185 

The  last  three  examples,  as  well  as  many  more  not  reported  here,  support 
the  theoretical  development  in  that  the  peak  estimator  <7^  has  a  smaller  MSE 
than  the  average  sum  of  squares  estimator  cr^  when  the  observations,  process 
values  and  peak  values,  constitute  random  samples.  The  tendency  to  have 
a  smaller  MSE  also  has  been  observed  more  generally  in  time  series  records 
when  the  sample  consists  of  the  entire  time  series.  This  is  illustrated  next. 
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4.2.2  Non-Random  Samples:  Discrete  Spectra 

In  the  preceding  subsection  we  dealt  with  random  samples  from  a  continuous 
spectrum  AR{2)  process.  In  this  subsection  we  continue  to  compare  and 
cr^  given  sinusoidal  data  plus  noise  of  the  form  (  8)  which  is  more  akin  to 
real  acoustic  emission. 

Our  sinusoidal  model  consists  of  4  Gaussian  sinusoids  plus  independent 
Gaussian  noise  of  the  form  (  8)  where  the  signal  to  noise  ratio  (SNR)  is 
IQdB.  We  have  generated  100  independent  time  series  from  the  process  each 
of  length  1000,  and  evaluated  <7^  and  <7^  from  each  series.  Given  a  time  series, 
0-2  is  evaluated  from  the  entire  record  (n  =  1000),  while  the  peak  estimator 

A 

(7^  is  obtained  from  all  the  peaks  in  the  series  greater  than  L.  This  gives  100 
cr^’s  and  100  (7^’s  from  which  we  obtained  an  average  and  a  sample  MSE.  The 
entire  procedure  was  repeated  12  times.  In  each  case  we  have  also  computed 
the  average  number  of  peaks  used  in  the  evaluation  of  the  cr^’s.  Observe  that 
on  the  average,  the  number  of  qualified  peaks  exceeding  L  is  much  smaller 
than  1000. 

As  can  be  seen  from  Tables  7  and  8,  in  this  particular  simulation  the 
peak  estimator  tends  to  underestimate  the  true  power,  hut  still,  as  in  random 
samples,  it  consistently  gives  smaller  MSE.  The  theory  of  statistics  judges 
the  overall  performance  of  estimators  by  the  smallness  of  their  MSE  (e.g.  see 
Bickel  and  Doksum  (1977),  p.  117). 
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Table  7.  Each  entry  gives  the  average  of  and  cr^  and  the  corresponding 
MSE  from  100  independent  sinusoidal  time  series.  L  =  0.9,  {u}i,U2,a}z,uj4)  = 


1.0, 1.1),  SNR 

=  10JB. 

The  true 

power  is  cr^ 

=  4.3969. 

No. 

a 

'■2 

A 

0-2 

Ave 

MSE 

Ave 

MSE 

Ave.  Peaks 

1 

4.42298 

4.16290 

4.01190 

3.72059 

148.010 

2 

4.34284 

4.53365 

3.91621 

3.98727 

148.200 

3 

4.41098 

3.68609 

3.96479 

3.37266 

148.170 

4 

4.50014 

3.55210 

4.01090 

3.04897 

151.350 

5 

4.35831 

4.37212 

3.91347 

3.85721 

148.030 

6 

4.49152 

4.39982 

3.99070 

3.88936 

151.040 

7 

4.33635 

3.63446 

3.86039 

3.36573 

150.130 

8 

4.69967 

5.38226 

4.20014 

4.22441 

150.170 

9 

4.44479 

3.75730 

4.04271 

3.36434 

147.660 

10 

4.46937 

4.54787 

4.00060 

4.00205 

149.130 

11 

4.39735 

4.84906 

3.97789 

4.49687 

149.030 

12 

4.53469 

3.96519 

4.10351 

3.54251 

149.160 

^  A 

Table  8.  Each  entry  gives  the  average  ofcr^  and  (7^  and  the  corresponding 
MSE  from  100  independent  sinusoidal  time  series.  L  =  1.3,  = 

(0.37,0.93,1.77,2.53),  SNR  =  IQdB.  The  true  power  is  cr'^  =  4.3969. 


No. 

(7 

2 

0-2 

Ave 

MSE 

Ave 

MSE 

Ave.  Peaks 

1 

4.17966 

3.90781 

3.26522 

3.75073 

192.050 

2 

4.53908 

4.34813 

3.55319 

3.38525 

195.600 

3 

4.61030 

3.70160 

3.63241 

3.01491 

197.230 

4 

4.06742 

4.15672 

3.22769 

4.04330 

184.870 

5 

4.50532 

2.99449 

3.53389 

2.69128 

198.240 

6 

4.57901 

4.35255 

3.63311 

3.53046 

191.790 

7 

4.51478 

4.38266 

3.57757 

3.65392 

189.690 

8 

4.50253 

4.17192 

3.59954 

3.43382 

187.180 

9 

4.59567 

4.81051 

3.65653 

3.52461 

187.800 

10 

4.91599 

6.33309 

3.86769 

4.13517 

194.010 

11 

4.57212 

4.69809 

3.62599 

4.02183 

191.800 

12 

4.62566 

5.45805 

3.65549 

4.09232 

189.130 
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5  Summary 

As  discussed  in  the  beginning  of  the  paper,  the  estimation  of  power  from 
noisy  acoustic  emission  time  series  is  important  for  source  characterization. 
In  this  regard,  we  have  illustrated-by  deriving  under  assumptions  a  class  of 
conditional  peak  distributions-the  fact  that  relatively  large  peaks  observed 
in  stationaxy  time  series  contain  information  about  the  power.  This  led  to 
the  construction  of  the  peak  estimator  (  21). 

The  power  can  be  estimated  straightforwardly  by  the  average  sum  of 
squares-the  sample  variance  when  the  mean  is  zero-or,  alternatively,  from 
the  peak  estimator  (21)  which  uses  peaks  exceeding  a  relatively  large  positive 
threshold  L.  Under  certain  conditions,  the  peak  estimator  tends  to  have  a 
smaller  MSE  than  the  average  sum  of  squares.  Simulations  show  that  the 
peak  estimator  may  underestimate  the  true  power  when  the  peaks  do  not 
form  a  random  sample  and/or  when  the  peaks  in  the  sample  are  not  large 
enough.  On  the  other  hand,  overestimation  may  occur  too  when  the  sample 
contains  extreme  peaks  and  L  is  not  large  enough. 

Common  sense  suggests  that  both  estimators  be  evaluated  and  used  in 
some  combination.  Since  under  some  conditions  spelled  out  in  the  Appendix 
the  average  sum  of  squares  is  erratic,  the  peak  estimator  is  a  quantity  to  be 
reckoned  with. 

When  troughs  occur  above  level  0  with  a  high  probability,  the  process  is 
not  narrow  band  to  a  sufficient  degree  and  L  must  be  relatively  large.  On  the 
other  hand,  an  excessively  large  L  will  clearly  disqualify  most  peaks  thus  forc¬ 
ing  the  evaluation  of  the  peak  estimator  from  unduly  small  samples.  This 
may  result  in  underestimation-when  the  sample  contains  too  many  values 
close  to  L-ox  in  overestimation-when  the  sample  contains  too  many  extreme 
peaks.  The  important  question  of  an  optimal  threshold  is  open  at  present. 
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6  Appendix:  Discrete  Spectrum  Considera¬ 
tion 

To  clarify  the  power  estimation  problem  in  mixed  and  discrete  spectrum 
processes,  we  consider  a  purely  discrete  spectrum  process.  For  this  we  follow 
Houdre  and  Kedem  (1994). 

Let  t  =  0,  ±1,  ±2,  •  •  •,  be  a  weakly  stationary  discrete  spectrum  pro¬ 

cess  with  mean  E[Xt]  =  m,  autocovariance  R(k)  =  E[(Xt+fc  —  m){Xt  —  m)], 
A:  =  0,  ±1,  ±2,  •  •  •,  and  spectral  distribution  F{du!)  supported  at  p-1-1  distinct 
atoms  in  {a;o,a;i,  •  •  •  ,uip]  €  (— 7r,7r],  such  that  ujq  =  0.  Then, 

Xt-m=  r 
J—K 

where, 

i=o  j=0 

and  the 

^(H})  = 

are  orthogonal  with  mean  zero.  Clearly  ^  :  B(— 7r,7r]  — v  L‘^{P)  is  <T-additive. 

As  an  estimator  for  R{k)  we  choose  the  (simplified)  sample  autocovari¬ 
ance 

4  -  m){Xi  -  m) 

<=i 

The  power  is  estimated  from  the  average  sum  of  squares 

1  ^ 

If  we  define 

li=0  J 

Then 

/27r 
■2ir 


25 


and  we  can  see  that  for  each  fixed  k,  the  lag  process  {Xi+k  —  fn){Xi  —  m) 
admits  a  Fourier  representation  with  respect  to  :  B(— 27r,27r]  — )■ 
and  a  dc  component 

>=0 

Theorem  6.1  Let 

(25)  =  m  +  <  =  0,±1,±2,--- 

j=o 

where  =  0,  and  for  j  ^  I  =  0,  6e  a  complex  valued  weakly  sta¬ 

tionary  process  with  mean  m,  autocovariance  R{-),  and  a  spectrum  supported 
at  p  \  distinct  atoms  {wq  =  0,  Wi,  •  •  •  ,LOp}  G  (— tt,  tt].  Then  as  N  —*  oo, 

(2«)  ^  -  ’n){X,  -  m)  i  =  0,  ±1,  ±2, . . . 

■'''  t=l 

Furthermore, 

(27)  ^  E{Xt+k  ~  m)(Xt  -  m)  R{k),  =  0,  ±1,  ±2,  •  •  • 

if  and  only  if 

(28)  =  ;  =  0,l,---,p 

where  the  <f)j  are  random  phases. 

The  requirement  of  constant  amplitudes  |^j|  means  that  if  {Xt)  in  (  25)  is 
a  real-valued  Gaussian  process,  its  sample  autocovariance  is  not  consistent. 
For  a  proof  of  Theorem  6.1  and  further  extensions  see  Houdre  and  Kedem 
(1994). 
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TIME  SERIES  FROM  "PLATE" 


Figure  1:  An  electrical  signal  generated  by  acoustic  emission  resulting  from 
ball  dropping  on  a  metal  plate. 


TIME  SERIES  FROM  'PLATE*  BLOCK' 


Figure  2:  An  electrical  signal  generated  by  acoustic  emission  resulting  from 
ball  dropping  on  a  metal  plate  plus  block. 


AVERAGED  FFT  FROM  THE  “PLATE*  REALIZATIONS 


Figure  3:  A  portion  of  average  FFT  from  9  independent  “plate”  realizations. 
n  =  4096,  At  =  200  x  10”^ 


averaged  FFT  FROM  THE  "PLATE  +  BLOCK"  REALIZATIONS 


Figure  4:  A  portion  of  average  FFT  from  10  independent  “plate+block”  re¬ 
alizations.  n  =  4096,  At  =  200  x  10“®. 


AVERAGE  AUTOCORRELATION  FOR  'PLATE'  REALIZATIONS 


Figure  5:  Average  autocorrelation  from  9  independent  “plate”  realizations. 
n  =  4096,  At  =  200  x  10"^ 


AVERAGE  AUTOCORRELATION  FOR  'PLATE+BLOCK'  REALIZATIONS 


Figure  6:  Average  autocorrelation  from  10  independent  “plate+block”  real¬ 
izations.  n  =  4096,  At  =  200  x  10“®. 


Acoustic  Emission  Energy 
(arbitrary  units) 


Figure  7:  Relationship  between  source  and  acoustic  emission  energies  for  a 
steel  ball  bearing  dropped  from  various  heights  onto  and  aluminum  alloy  plate. 
The  sample  correlation  coefficient  is  0.974- 


